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1.  Introduction 


The  original  goal  of  radar  (radio  detection  and  ranging)  was  to  provide  informa¬ 
tion  about  the  existence  and  position  of  targets.  It  was  soon  realized  that  additional 
capabilities  were  desirable,  in  particular  the  ability  to  classify  the  target.  Depending 
on  the  nature  of  the  application,  classification  might  be  concerned  with  friend  or  foe 
status  or  with  flight  rule  status  ( visual /instrumental)  or  many  other  things.  In  some 
cases  it  is  possible  to  obtain  this  information  by  requiring  the  aircraft  to  transmit 
it  via  transponders  or  IFF  devices;  in  other  cases  this  is  impossible  or  undesirable. 
The  latter  problem  is  the  so  called  non-cooperative  radar  target  identification  (RTI) 
problem. 

From  physical  considerations  it  is  known  that  information  about  the  shape  of 
the  target  is  present  in  the  radar  return  signal,  since  the  target  will  distort  the 
incoming  radar  signal  in  a  geometry-dependent  manner.  The  problem  of  RTI  is  to 
extract  useful  information  from  the  distortions  in  the  radar  return.  This  extraction 
process  can  be  considered  to  consist  of  two  stages:  first  there  is  a  signal  processing 
stage,  where  the  radar  return  signal  is  processed  to  extract  some  discriminating 
information;  second  is  the  decision  stage,  where  the  extracted  information  from 
the  signal  processing  stage  is  matched  against  a  data  base  of  known  targets  (the 
catalog),  and  a  decision  as  to  the  nature  or  identification  of  the  target  is  made. 

One  method  of  performing  target  identification  is  to  use  no  signal  processing 
step,  so  the  decision  algorithm  operates  directly  on  the  frequency  data.  This  method 
has  the  advantages  of  simpler  implementation  and  of  not  losing  information  in  the 
signal  processing  stage.  However,  there  are  some  potential  disadvantages.  First,  if 
the  number  of  frequency  measurements  is  large,  the  decision  step  will  be  computa- 
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tionally  intense.  Second,  the  frequency  domain  information  does  not  relate  to  the 
geometric  configuration  of  the  target  in  a  direct  manner;  thus,  one  cannot  directly 
obtain  geometric  information,  such  as  overall  length,  about  the  target. 

This  report  considers  an  alternative  approach  to  target  identification.  In  this 
approach,  the  signal  processing  step  is  designed  to  extract  geometrically  relevant 
scattering  characteristics  of  the  target.  This  method  of  RTI  consists  of  a  signal 
processing  step  which  converts  the  radar  return  from  the  frequency  to  the  time 
domain,  to  obtain  the  radar  cross-section  versus  time  signal  of  the  target.  The 
advantage  of  this  method  is  that  the  relation  with  the  geometry  of  the  target  is 
much  more  direct  in  this  representation.  The  decision  step  then  classifies  the  target 
from  the  time  domain  information.  Traditionally  the  transition  to  the  time  domain 
has  been  accomplished  using  Fast  Fourier  Transform  (FFT)  techniques  [1].  FFT 
techniques  suffer  from  two  potential  disadvantages,  however.  First,  they  do  not  yield 
any  data  reduction,  so  the  output  of  the  signal  processing  stage  gives  at  least  as 
many  data  points  in  the  time  domain  as  originally  present  in  the  frequency  domain. 
Thus,  there  is  no  decrease  in  processing  requirements  for  the  decision  algorithm. 
Second,  the  time  resolution  of  the  FFT  is  limited  by  the  bandwidth  of  the  frequency 
domain  return,  so  closely  spaced  features  on  the  target  might  not  be  resolved. 

This  report  addresses  the  option  of  using  modern  signal  processing  methods 
to  convert  from  the  frequency  to  the  time  domain.  In  particular.  Autoregressive 
Moving  Average  (ARMA)  models  are  used  to  parameterize  the  radar  data.  ARMA 
models  do  not  suffer  from  the  resolution  limit  of  the  FFT  methods.  In  addition, 
ARMA  models  are  parametric,  so  the  information  present  in  a  large  number  of 
frequency  samples  can  be  concentrated  into  a  small  number  of  parameters,  thus 
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reducing  the  computational  burden  for  the  decision  stage.  Finally,  the  parameters 
in  the  ARMA  model  axe  closely  related  to  relevant  geometric  features  of  the  target. 
As  a  result,  even  if  the  target  cannot  be  identified,  some  geometric  information 
about  the  target  can  be  extracted  from  the  ARMA  parameters. 

In  this  report,  we  study  the  application  of  ARMA  models  to  the  signal  processing 
part  of  the  RTI  problem.  First,  we  develop  a  special  representation  for  the  ARMA 
model  suited  to  the  specifics  of  the  RTI  problem.  We  then  describe  an  algorithm  for 
estimating  the  needed  parameters  from  the  frequency  data.  We  study  the  effects  of 
several  improvements  over  the  basic  algorithm,  including  the  use  of  singular  value 
decomposition  and  the  elimination  of  spurious  scattering  centers  from  the  estimated 
model.  Both  lead  to  an  algorithm  that  is  more  robust.  Finally,  we  present  results  of 
this  algorithm  are  shown  using  radar  data  obtained  by  compact  range  measurements 
of  commercial  aircraft  models. 
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2.  Canonical  Scatterers  and  ARMA  Models 


This  chapter  discusses  the  modeling  of  radar  targets  by  canonical  scattering 
centers,  and  the  relation  to  ARMA  models.  In  particular,  it  will  be  shown  that 
ARMA  models  are  an  appropriate  way  of  modeling  radar  returns. 

2.1  Canonical  Scatterers 

It  is  well-known  [2]  that  scattering  of  electromagnetic  waves  at  high  frequencies, 
i.e.  at  frequencies  such  that  the  wavelength  is  small  compared  to  the  scatterer,  can 
be  modeled  as  a  set  of  independent  scattering  centers.  There  are  several  varieties 
of  canonical  scattering  centers  used  in  this  process;  typical  examples  are  [2]  the 
point,  edge  and  plate  scatterers  and  the  comer  reflector.  As  is  stated  in  [2],  these 
examples  are  typical  of  the  types  of  scatters  that  can  be  found  on  radar  targets. 

From  [3],  we  have  that  all  the  above  canonical  types  of  scatterers  can  be  modeled 
in  the  following  way:  if  a  target  consisting  of  m  scatterers  is  illuminated  using  a 
CW  radar,  then  the  backscattered  field  received  by  a  linearly  polarized  antenna  will 
give  a  signal  of  the  form: 

yk  =  t,Kiez‘^W*h)'-  (2.1) 

1=1 

where  K{  is  the  amplitude,  r,  the  relative  range,  and  #,■  the  type  of  the  scatterer; 
fk  is  the  k™  frequency,  and  m  is  the  number  of  scatterers.  For  the  canonical 
scatterers  in  [3],  t,  takes  on  the  following  values 

ti  G  {  — 1, —0.5, 0,0.5, 1}  (2.2) 
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In  addition,  for  stepped  frequencies,  we  can  represent  fk  by 

fk  =  fo  +  k6f,  k  =  0,l,2,...,N-l  (2.3) 

where  /0,  8/  and  N  are  given  parameters  of  the  system  under  consideration.  The 
frequency  /0  represents  the  lowest  frequency  the  radar  system  transmits  and  8j 
is  the  size  of  the  frequency  step.  Estimation  of  the  parann  n  is  in  (2.1)  from  N 
measurements  leads  to  a  canonical  scattering  description  of  the  target;  see  [3].  Un¬ 
fortunately,  this  estimation  is  nonlinear,  and  is  not  practical  for  real-time  target 
identification  at  this  time. 

For  reasonably  small  bandwidths,  a  feasible  method  of  target  identification  can 
be  obtained  by  approximating  the  canonical  scattering  model  of  (2.1).  In  particular, 
the  factor  (j2nfk)u  can  be  approximated  by  a  factor  of  the  form  pf  for  some  pt. 
This  approximation  yields  the  data  model: 

m 

n  =  (24) 

1=1 

where  d{  is  in  general  not  equal  to  I\,  in  (2.1)  because  of  the  approximation  involved. 
The  variable  p;  can  be  found  from: 

2ri2wfl 

Pi  =  3  c  (2.5) 

or,  conversely,  we  can  find  pi  and  r,  from  p  bv: 


Pi  =  |Pi| 

*  “  <5-^ 


(2.6) 

(2.7) 


where  R  =  is  the  maximum  unambiguous  range  of  the  radar.  In  practice  the 
absolute  range  is  not  used;  instead,  r,  represents  the  range  relative  to  a  zero  range 


reference  point. 

The  physical  relevance  of  the  pi  parameter  is  dubious  because  it  is  defined  from 
an  approximation.  However,  it  is  important  to  estimate  />,  since  it  allows  the  calcu¬ 
lation  of  a  physically  important  quantity:  the  amount  of  energy  associated  with  a 
scattering  center.  This  energy  can  found  by  considering  the  contribution  of  a  single 
scattering  center  to  the  total  signal.  The  single  scattering  center  component  is 

{diPi}k=Q 

The  energy  is  the  square  of  the  absolute  value  of  this  contribution,  or 

N-l  i  _  2 N 

P,  =  E  IAf  M"  =  (2-8) 

fc=0  1  Pi 

Equations  (2.4)  and  (2.8)  will  be  used  to  form  a  parametric  model  of  target  scat¬ 
tering. 

2.2  The  ARMA  Modeling  Interpretation 

In  the  previous  section  we  established  that  the  return  of  a  stepped  frequency 
radar  can  be  modeled  as 

m 

yk  =  YldiPi 

t=i 

In  this  section  we  will  show  how  this  relates  to  Autoregressive  Moving  Average 
(ARMA)  modeling  and  the  radar  cross  section  versus  time  relation  of  the  target. 

If  we  take  the  inverse  Fourier  Transform  of  the  above  relation  for  y*,  v>v  .  ^t 

m  d 

yw of-m  (2-9> 
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where  r  is  relative  range.  The  square  of  the  absolute  value  of  this  function  gives 
the  radar  cross  section  as  a  function  of  relative  range. 

If  we  change  from  the  Fourier  to  2  Transform  (i.e.  z  —  eJ7r(1_2r/ft) );  equation 
(2.9)  becomes 


no 


dj 

Pi 


C(z) 

A{z) 


(2.10) 


where 

N  A(z)  —  IK^r  -p^  =  zp  +  ai zp~x  +  ...  +  dP  (2-11) 

i=i 

and  where  C(z)  is  a  polynomial  of  degree  less  than  or  equal  to  m  —  1.  It  can  be 
seen  that  the  p,  coefficients  axe  the  zeros  of  A(z),  and  so  they  are  the  the  poles 
of  Y(z).  Hence  the  model  we  have  assumed  for  the  radar  data  corresponds  to  an 
autoregressive  moving  average  model  with  m  poles  and  m  —  1  zeros. 

Note  that  the  exponential  model  (2.4)  applies  to  the  frequency  domain,  and 
the  rational  model  (2.9),  (2.10)  applies  to  the  time  domain.  This  is  in  contrast  to 
standard  time  series  analysis  treatments  of  ARMA  models,  where  the  damped  ex¬ 
ponential  model  applies  to  the  time  domain,  and  the  corresponding  rational  transfer 
function  model  is  in  the  frequency  domain  [4].  This  difference  is  cosmetic  only,  and 
is  easily  handled. 

Comparing  (2.4)  with  (2.10),  it  can  be  seen  that  each  term  {d-pf}^1  corre¬ 
sponds  to  a  term  As  has  been  stated  above,  the  terms  {dj?*]  N  '}  corn  pond 
to  the  scattering  centers  of  the  target;  hence,  the  terms  model  the  r<  p-nse 
of  the  individual  scattering  centers  in  the  time  domain.  In  other  words,  i  !■•  'ime 
domain  model  (2.10)  describes  the  radar  cross  section  (RCS)  versus  time  response 
of  the  target  as  the  sum  of  the  RCS  responses  of  the  scattering  centers.  The  RCS 
time  response  can  be  obtained  by  evaluating  (2.9)  for  0  <  r  <  R  or  equivalently, 
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evaluating  (2.10)  for  z  =  e^,  0  <  <f>  <  2n.  This  RCS  time  response  is  similar  to  the 
responses  obtained  by  taking  the  DFT  of  the  frequency  data,  but  does  not  have  the 
resolution  limitation  inherent  in  the  DFT. 

The  nature  of  each  scattering  center  is  determined  by  the  values  of  the  param¬ 
eters  pi  =  p,eJir(1-2r*^  and  d,  associated  with  it.  First,  rx  determines  the  relative 
range.  Second  the  magnitude  p,  of  p,  determines  the  shape  of  the  RCS  response 
curve;  if  pi  is  close  to  one  (i.e.  p<  is  close  to  the  unit  circle),  the  curve  will  be  a  nar¬ 
row  and  sharp  peak,  indicating  that  the  scattering  is  concentrated  at  a  particular 
range.  If  on  the  other  hand,  p,-  is  not  close  to  one  (i.e.  p,  is  far  away  from  the  unit 
circle),  the  curve  is  much  wider  and  has  no  sharp  peak,  indicating  the  scattering 
center  response  is  spread  out  in  range.  Finally,  the  values  of  pi  and  di  combine  to 
give  the  energy  P,  associated  with  the  scattering  center,  as  indicated  in  (2.8). 

The  above  ideas  are  illustrated  in  Figure  2.1.  This  Figure  shows  the  response 
|F(r)|  of  a  hypothetical  target  composed  of  four  scattering  centers.  The  lower 
graph  shows  the  location  of  the  zeros  {p,}f_:1  in  the  complex  plane.  By  comparing 
the  locations  of  the  zeros  and  the  positions  of  the  peaks  in  the  response  graph,  we 
can  clearly  see  how  the  argument  of  p,  relates  to  the  range  of  the  corresponding 
peak.  Also,  by  comparing  peaks  3  and  4,  we  can  see  that  when  p,  is  closer  to  the 
unit  circle,  the  corresponding  peak  will  be  sharper  and  narrower.  Peaks  1  and  2 
show  that  even  if  two  scattering  centers  have  pi’s  that  are  equally  close  to  the  unit 
circle,  the  amplitude  of  their  peaks  may  vary  greatly,  since  the  d,  values  may  be 
different.  For  reference,  Figure  2.2  shows  the  FFT  of  the  same  hypothetical  target. 
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3.  Techniques  to  Estimate  the  Model  Parameters  from  Radar  Data 


Based  on  the  results  of  the  previous  section,  the  frequency  domain  radar  backseat  - 
ter  data  can  be  modeled  as  a  sum  of  complex  exponential  terms.  More  practically, 
the  data  will  also  contain  an  additive  noise  component  due  to  receiver  noise,  clut¬ 
ter,  etc.  In  either  case,  the  signal  processing  objective  is  to  estimate  the  parameters 
which  describe  the  exponential  components  from  a  finite  set  of  (possibly  noisy) 
data.  This  problem  is  considered  in  the  time  series  analysis  literature  as  the  so- 
called  Prony  estimation  problem  [4]. 

Below,  we  develop  algorithms  which  estimate  the  needed  parameters  from  radar 
data.  These  algorithms  are  based  on  singular  value  decomposition  (SVD)  versions 
of  Prony  algorithms.  We  then  discuss  modifications  to  the  algorithms  to  make  tin  m 
more  appropriate  for  radar  data. 

3.1  The  Backward  Prediction  Equations 

Assume  that  we  wish  to  model  the  radar  data  as  a  sum  of  complex  exponentials, 
as  in  equation  (2.4),  and  that  the  measured  data  contains  noise:  thus 

m 

yk  =  YldiPi  +  k  =  0, . . .  ,N  —  l  (3.1) 

i=i 

where  {n^  is  a  sequence  of  zero  mean  noise.  For  concreteness,  we  will  assume 
Ufc  is  white  noise  with  zero  mean  and  variance  A2;  however,  the  colored  noise  case 
can  also  be  handled,  using  algorithms  similar  to  the  ones  we  describe  below  [5]. 
Below  we  discuss  various  alternative  ways  to  estimate  the  px  and  dt  coefficients 
from  the  measured  data  sequence  {yk}k=o  ■  These  methods  are  adapted  from  similar 
algorithms  discussed  in  [4]. 
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First  we  consider  the  estimation  of  the  p,  parameters.  Recall  that  the  p,  param¬ 
eters  are  the  zeros  of  the  A(z )  polynomial  (see  equation  (2.10)).  It  can  be  shown 
[4]  that  the  coefficients  a,  of  the  A(z)  polynomial  satisfy  the  equation: 

m 

yk  +  =  ° 

i=l 


This  can  be  written  in  matrix  form  as: 


Vm+i  Vm 
Vm+2  Vm+i 

Vn-  1  VN-2 


(3.2) 


Where  M  is  some  integer  at  least  equal  to  m  and  where  equality  holds  in  the 
noiseless  case;  we  will  discuss  the  selection  of  M  below.  Solutions  based  on  (3.2) 
are  called  total  least  squares  solutions  [6j.  It  was  shown  in  [6]  that  total  least 
squares  solutions  generally  have  superior  accuracy  when  compared  with  ordinary 
least  squares  methods  [6]. 

It  should  be  noted  here  that  in  a  given  application  we  would  not  know  the  correct 
value  of  m.  This  point  will  be  further  discussed  below. 


3.2  Use  of  the  Singular  Value  Decomposition 

Estimation  of  parameters  from  a  short  data  record  is  a  difficult  problem;  while 
most  popular  techniques  are  quite  good  at  estimating  parameters  from  noisy  data 
if  the  record  length  is  large,  their  ability  to  perform  well  for  short  records  is  signif¬ 
icantly  less.  One  technique  modem  signal  processing  research  has  found  to  combat 
this  effect  is  the  use  of  the  Singular  Value  Decomposition  (SVD).  The  SVD  is  an 
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orthogonal  decomposition  technique  for  matrices;  it  can  be  regarded  as  a  general¬ 
ization  of  the  eigenvalue  decompositions  that  exist  for  square  matrices. 

The  use  of  the  SVD  to  the  estimation  problem  lies  in  its  ability  to  help  determine 
the  rank  of  a  noisy  matrix.  Recall  that  estimating  the  polynomial  A(z )  requires 
forming  a  matrix  out  of  the  noisy  data  (equation  (3.2)).  In  the  noiseless  case,  the 
rank  of  this  matrix  is  m,  the  number  of  scattering  centers  on  the  target.  However, 
due  to  the  effect  of  the  noise  on  the  data,  the  actual  rank  will  be  greater  than  m. 
For  small  noise  variances,  the  data  matrix  will  be  “close”  to  a  matrix  of  rank  m,  and 
the  SVD  can  help  find  the  closest  rank  m  approximation.  It  is  expected  that  this 
procedure  will  reduce  the  effects  of  the  noise  on  the  estimate.  For  a  more  detailed 
discussion,  see  [7,8,9]. 

Next  we  need  to  consider  how  to  implement  the  SVD  technique  in  our  algorithm. 
To  do  this  we  modify  equation  (3.2)  as  follows:  we  first  compute  the  matrix  on  the 
left  of  (3.2)  and  then  compute  its  SVD: 

Y  =  UY.VH  (3.3) 

where  Y  is  the  data  matrix,  U  and  V  are  unitary  matrices,  and  E  is  a  diagonal 
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matrix,  given  by: 


<7  X  0  ...  0 

0  <72  .  .  .  0 

0  0  ...  <tm+i 

00...  0 


(3.4) 


0  0 


0 


J  (N-M-l)x(M+l) 


where 


cri  >  02  >  •  •  •  >  1  >  o.  (3.5) 

As  noted  above,  if  the  SNR  is  sufficiently  high,  we  will  have: 


0\  —  •  •  •  —  ^  *^m+ 1  —  •  •  •  ^ 


(3.6) 


where  <7m+1 . . .  crM+1  are  nonzero  only  because  of  the  nonzero  noise  term.  Thus, 
some  effects  of  the  noise  can  be  eliminated  by  replacing  S  by  S'  defined  as 


S'  = 


°\ 


0 


<*n 

0  0 


J  (N-Af-1)x(M+l) 


(3.7) 


for  some  value  of  n  satisfying  m  <  n  <  M  +  1  (we  will  discuss  ways  of  choosing  n 
below).  Finally,  the  vector  a  =  [ldi . . .  a\f]T  is  found  by  solving 


UZ'VHaM  =  0 


(3.8) 
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in  the  least  squares  sense  [6]. 

The  SVD  also  greatly  helps  the  order  estimation  problem.  Assume  we  do  not 
know  the  true  number  of  scattering  centers  m,  but  are  confident  that  the  value  of 
M  is  large  enough  to  model  the  scattering  behaviour  of  the  target.  Then  we  can 
select  n  <  M  to  best  model  the  target  and  reject  noise.  Given  the  short  data  lengths 
typical  of  RTI  problems,  a  large  enough  M  can  typically  be  achieved  by  choosing  M 
to  be  a  large  fraction  of,  or  even  equal  to,  y.  Then  n  can  be  chosen  by  considering 
the  distribution  of  the  sv’s:  we  choose  n  such  that  <7n+\  is  the  first  sv  that  falls 
below  a  preset  threshold.  In  this  way  <j\ . . .  an  represent  the  important  (high  energy) 
signal  components,  and  the  other  sv’s  represent  either  small  signal  components  or 
noise.  Note  that  we  may  not  model  all  of  the  target  scattering  components,  as 
the  weaker  ones  can  be  eliminated  from  the  model  by  choosing  n  small;  on  the 
other  hand,  the  weaker  scattering  components  are  more  prone  to  be  misestimated 
in  noise.  Thus,  this  procedure  represents  a  trade-off  between  geometric  faithfulness 
and  robustness  in  the  presence  of  noise.  The  clear  advantage  of  this  method  is 
its  built-in  robustness:  it  selects  the  high  energy  signal  components,  which  are 
expected  to  suffer  the  least  degradation  from  the  noise.  Another  advantage  of 
this  method  over  order  selection  in  a  straightforward,  non-SVD  estimator  is  that 
because  the  order  reduction  proceeds  in  an  orthogonal  manner,  we  do  indeed  select 
the  strongest  signal  components,  whereas  the  straightforward  algorithm  will  select 
mixtures  of  several  components,  as  can  be  seen  from  the  experiments  in  [10].  The 
SVD  method  is  one  of  the  best  currently  known  practical  methods  for  tackling  the 
order  estimation  problem  for  radar  data. 
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3.3  Estimation  of  Modal  Amplitudes 


Once  we  have  estimates  of  the  A(z )  polynomial  coefficients,  we  can  find  the 
corresponding  zeros  p,  by  using  a  polynomial  root  finding  method.  It  then  remains 
to  find  the  corresponding  amplitude  ( d{ )  coefficients.  From  equation  (2.4)  we  can 
see  that  given  the  p,,  the  problem  of  estimating  the  d,  is  a  linear  one,  so  we  can  use 
a  standard  Least  Squares  approach.  Of  course  we  do  not  know  the  correct  values 
of  the  pi,  but  we  can  use  the  estimates  obtained  as  described  above  in  their  stead. 
The  least  squares  approach  to  estimating  the  d,  coefficients  can  be  formulated  as 
follows: 


'  f?  •• 

rP 

Pm 

d0 

Vo 

.. 

~N-1 

Pm  \ 

1 

VN- 1 

(3.9) 


where  are  the  estimates  obtained  from  the  first  step.  This  equation  is  simply 
(3.1)  written  out  in  matrix  form.  Since  the  equation  (3.9)  is  overdetermined,  we 
can  solve  it  in  a  least  squares  sense. 


3.4  Elimination  of  Spurious  Zeros 

In  previous  sections  we  computed  an  estimate  of  a  polynomial  Am(z )  of  order 
M  >  m,  and  a  corresponding  set  of  M  zeros  {p,}^.  We  also  computed  ampli¬ 
tudes  corresponding  to  these  zeros.  The  next  step  is  to  separate  the  m  zeros  that 
correspond  to  the  radar  target  (signal  zeros)  from  the  remaining  spurious  zeros. 

In  [10],  we  reported  a  criterion  for  zero  inspection  that  worked  well  with  a 
related  ARMA  algorithm.  This  criterion  is  based  on  the  observation  that  a  zero 
with  modulus  greatly  different  from  one  will  influence  only  the  data  elements  at  the 
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start  (if  the  modulus  is  smaller  than  one)  or  at  the  end  (if  it  is  greater  than  one). 
In  the  time  domain,  we  observe  that  such  zeros  correspond  to  widely  distributed 
scattering  centers,  which  are  of  little  interest  to  radar  target  identification.  As 
reported  in  [10],  removal  of  such  zeros  improved  the  performance  of  the  algorithm 
we  were  then  using. 

Unfortunately,  the  criterion  in  [10]  does  not  seem  to  mesh  well  with  the  current 
SVD-based  algorithm.  Results  of  applying  this  zero  selection  method  on  aircraft 
data  shows  two  undesirable  characteristics.  First,  this  criterion  produces  erratic 
changes  in  the  selected  zeros  as  one  changes  the  value  of  m.  This  is  undesirable 
for  RTI  applications  because  our  experience  has  shown  that  one  needs  to  select  the 
order  differently  for  different  aircraft.  Therefore,  it  is  necessary  for  the  order  to  be 
selected  on-line,  and  thus  we  need  algorithms  that  will  produce  similar  estimates 
even  if  the  order  selected  for  the  noisy  data  differs  from  the  order  selected  for  the 
catalog  prototype.  Second,  the  experimentally  found  spread  of  the  resulting  d,  and 
Pi  parameters  was  such  as  to  make  them  useless  for  classification:  there  were  no 
clear  clusters  for  the  different  targets.  Hence,  we  needed  a  better  way  to  differentiate 
spurious  zeros  from  the  signal  zeros.  We  are  currently  investigating  two  different 
techniques,  as  detailed  below. 

The  first  technique  is  a  refinement  of  the  method  in  [10].  Whereas  in  [10]  we 
kept  all  zeros  within  an  ring  around  the  unit  circle  with  outer  radius  r  and  inner 
radius  *  ,  we  now  choose  the  outer  and  inner  radii  of  the  ring  independently.  In 
particular,  since  one  expects  spurious  zeros  to  lie  inside  the  unit  circle  when  using 
backward  prediction  equations  [4],  the  inner  radius  is  set  much  closer  to  the  unit 
circle  than  the  outer  one. 
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The  second  technique  relies  on  the  observation  that  since  the  spurious  zeros 
correspond  to  modes  that  are  not  present  in  the  data,  their  estimated  energy  should 
be  low.  To  take  advantage  of  this  property,  the  algorithm  is  modified  as  follows. 
First,  all  M  zeros  are  kept,  and  the  amplitudes  {d,}^j,  corresponding  to  those  zeros 
are  found  using  equation  (3.9).  Then,  the  energy  of  each  mode  is  computed  using 
equation  (2.8).  The  m  zeros  whose  corresponding  energy  is  the  highest  are  kept 
and  the  remaining  ones  are  discarded.  Then,  the  energy  estimates  for  the  m  kept 
zeros  can  be  refined  by  recomputing  the  amplitude  estimation  using  (3.9)  with  only 
the  m  selected  zeros,  however,  since  the  energies  for  the  discarded  zeros  is  low,  the 
difference  between  the  two  estimates  should  generally  be  small,  so  this  refinement 
step  may  be  omitted. 

Preliminary  simulation  results  indicate  that  the  zero  selection  technique  which 
uses  the  energy  criterion  performs  significantly  better  than  the  one  using  zero  radii. 
This  will  be  seen  in  the  experiments  presented  in  section  4. 

Table  3.1  summarizes  the  ARMA  modeling  algorithm  used  in  this  report. 

3.5  Classification  Aspects  1 

In  this  section  we  report  results  pertaining  to  experiments  we  performed  to  get 
more  insight  into  the  ability  of  the  algorithms  discussed  to  form  the  signal  processing, 
stage  before  a  classifier.  We  have  looked  at  the  experimentally  found  distribution 
function  for  the  distance  of  the  two  strongest  scattering  centers  on  a  Boeing  747, 
at  the  clustering  properties  of  the  p,  and  parameters,  and  at  the  residual  error 
(defined  below)  as  a  function  of  n. 

^he  authors  would  like  to  thank  O  S.  Sands  for  his  contributions  to  this  material. 
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Table  3.1:  Summary  of  the  Algorithm. 


1.  Find  the  zeros  {p,}£fj  from  the  frequency  data  { yt }  £=0  as  fol¬ 
lows: 

(a)  Given  M,  form  the  matrix  Y  in  (3.2)  and  obtain  a  singular 
value  decomposition  of  Y  as  defined  (3.3). 

(b)  Given  a  cutoff  criterion,  find  the  smallest  n  such  that  ern 
satisfies  the  criterion,  but  crn+1  does  not.  Now  solve  the 
equation 

UY,'VHaM  =  0 

for  a.M ,  where  £7,  S'  and  V  axe  defined  in  (3.3),  (3.7)  and 
(3.8). 

(c)  Find  the  zeros  {pi}iix  of  the  polynomial  Am(z). 

2.  Find  the  modal  amplitudes  and  energies  {Pi}*i i  associ¬ 

ated  with  {pi)iLx  from  equations  (3.9)  and  (2.8) 

3.  Select  the  m  zeros  with  the  highest  energies  Pi.  From  these  m 
zeros  {pt}£Li  compute  the  ranges  {rJJ’lj  and  dispersions  {pi}'ll1 
using  equations  (2.7)  and  (2.6).  Optionally  recompute  the  m 
amplitudes  and  and  m  energies  from  equations  (3.9)  and  (2.8) 
using  the  m  selected  zeros. 

4.  The  target  response  Y(r)  can  be  found  from  (2.9)  using  the 

coefficients  {rt,  p<,  . 
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The  distribution  of  the  difference  in  relative  range  between  the  two  scattering 
centers  with  the  highest  energies  (i.e.  r\  —  r2)  for  the  Boeing  747  in  5  dBsm  noise  is 
clearly  Gaussian  to  good  approximation  in  the  case  of  the  algorithm  as  presented  in 
Table  3.1,  and  it  is  clearly  not  Gaussian  for  the  case  of  the  algorithm  with  the  ring 
criterion  described  in  the  previous  section  (see  Figure  3.1).  Moreover,  in  the  former 
case  the  variance  is  lower.  This  is  important  because  it  is  easier  to  design  classifiers 
for  parameters  with  a  Gaussian  distribution  than  for  for  other  parameters.  The 
most  important  point  however,  is  that  the  ring  criterion  algorithm  shows  significant 
bias,  whereas  the  algorithm  as  presented  in  Table  3.1  does  not.  The  algorithm  as 
presented  in  the  table  gives  a  mean  value  for  r2  —  rq  that  corresponds  to  the  relative 
range  found  in  the  noiseless  case;  as  can  be  seen  from  Figure  3.1,  the  ring  algorithm 
gives  a  mean  that  differs  from  the  true  mean  by  more  than  10%. 

It  has  further  been  found  that  the  values  of  the  dt  and  p,  coefficients  tend  to 
cluster  in  different  places  in  the  d  x  p  plane  for  the  different  aircraft;  this  means 
that  these  parameters  carry  information  that  will  aid  classification. 

Finally,  we  have  found  that  for  the  ring  criterion  algorithm  the  residual  <  rror  is 
not  a  monotonic  function  of  the  number  of  singular  values.  The  residual  error  is 
defined  as: 


N 


EK?  (v*  -  TZx  diPf)2 


(3.10) 


V'A'-l  „ 

2->fc= 0  Vk 

It  is  a  measure  of  how  well  the  estimate  models  the  data.  It  is  known  that  without 


the  elimination  of  spurious  zeros,  Ere^m)  will  decrease  monotonically  as  a  function 
of  the  number  of  singular  values  kept,  n.  This  is  an  important  property,  since 
it  guarantees  that  increasing  the  number  of  singular  values  can  only  improve  the 
quality  of  the  model.  It  was  found  experimentally  however,  that  the  ring  criterion 
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Relative  range  between  top  two  scatterers 


Figure  3.1:  Relative  Range  between  the  the  two  strongest  scatterers  on  a  Boeing 
747,  using  20  data  points  from  6.25-12  GHz,  10th  order  model,  7  sv’s  kept. 


algorithm  did  not  have  this  property.  In  fact,  for  this  algorithm  there  was  no 
predictable  regularity  to  the  error  versus  singular  values  function.  For  the  purposes 
of  RTI,  it  is  important  that  spurious  zeros  be  eliminated,  so  as  to  avoid  confusing 
the  classifier  with  physically  irrelevant  data.  It  was  for  this  reason  that  we  started 
to  investigate  the  energy  based  criterion  for  keeping  zeros,  and  indeed,  we  found 
that  for  the  algorithm  presented  in  Table  3.1,  the  residual  error  is  a  monotonically 
decreasing  function  of  the  number  of  singular  values  kept,  as  desired. 
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4.  Estimation  Results  Using  Compact  Range  Data 


This  section  describes  the  results  of  the  signal  processing  procedure  detailed  in 
the  previous  section  when  applied  to  compact  range  measurements  of  targets.  The 
data  used  contain  stepped  frequency  measurements  of  models  of  various  commercial 
aircraft;  specifically  the  DC  10,  the  Concorde  and  the  Boeings  707,  727  and  747 
are  compared.  The  measurements  have  all  been  made  at  the  OSU  ElectroScience 
Laboratory  Compact  Range,  and  are  available  in  various  frequency  ranges  from  2- 
36  GHz.  The  size  scales  for  the  model  aircraft  range  from  130-200,  so  the  frequencies 
correspond  to  the  5-280  MHz  for  a  full-size  aircraft,  depending  on  the  scale  of  the 
of  the  aircraft  under  consideration. 

The  stepped  frequency  measurements  consist  of  amplitude  and  phase  of  the 
radar  return  from  the  target  at  each  frequency.  The  radar  receiver  hardware  range 
gates  the  returning  signal  to  eliminate  spurious  scattering  effects  from  the  anechoic 
chamber.  The  data  is  then  calibrated  by  subtracting  the  return  of  the  empty  cham¬ 
ber  (with  pedestal)  and  dividing  by  the  calibrated  return  for  a  reference  target.  The 
calibrated  measurements  obtained  are  then  essentially  error-free  [11],  For  those  ex¬ 
periments  where  we  have  used  noisy  data,  we  have  added  white  Gaussian  noise  to 
these  measurements. 

4.1  Tests  of  Sensitivity  to  Parameter  Selection 

First  we  look  at  the  sensitivity  of  the  ARMA  algorithm  with  respect  to  choices 
of  certain  parameters  in  the  algorithm,  such  as  model  order  M,  number  of  singular 
values  kept  n,  etc.  We  also  consider  this  sensitivity  for  the  two  methods  of  zero 
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selection  as  described  in  Section  3.4.  We  first  show  the  algorithm  sensitive  to 
parameter  selections  when  a  radius  test  is  used  for  zero  selection. 

If  the  radius  test  method  is  used  (see  Section  3.4),  then  Step  3  in  Table  3.1  is  not 
performed,  and  instead  in  Step  lc  the  value  m  is  determined  as  the  number  of  zeros 
that  fall  within  a  preset  ring  of  inner  radius  rj  <  1  and  outer  radius  r2  >  1 .  These 
m  zeros  are  selected  and  all  following  steps  are  performed  with  only  this  subset  of 
zeros.  This  implies  that  all  subsequent  steps  use  m  instead  of  M . 

Figure  4.1  shows  the  estimated  response  for  the  Concorde  model,  using  model 
frequencies  from  18  to  26  GHz  (this  corresponds  to  frequencies  from  138-200  MHz 
on  the  full-size  Concorde),  N  =  30  data  points,  an  order  M  =  10,  while  keeping 
n  =  5  singular  values  and  m  =  2  zeros.  The  same  figure  also  shows  a  response 
curve  that  corresponds  to  the  same  situation  except  that  only  3  sv’s  are  kept  in 
the  estimate.  For  reference,  the  FFT  of  the  Concorde  data  is  shown  Figure  4.2. 
This  figure  shows  the  absolute  value  (in  dB)  of  the  inverse  Fourier  Transform  of 
the  data  used  by  the  algorithm.  The  inverse  Fourier  Transform  has  been  evaluated 
by  zero  padding  the  30  data  points  to  512  and  taking  the  IFFT.  It  is  clear  that 
the  ARMA  algorithm  is  sensitive  to  changes  in  the  number  of  singular  values  kept. 
Similar  observations  can  be  made  for  the  number  of  zeros  kept  and  the  model  order. 
This  verifies  the  claim  that  the  ring  method  of  zero  selection  is  sensitive  parameter 
choices.  The  figure  given  is  a  typical  case  taken  out  of  a  large  number  of  trials  that 
we  have  performed. 

If,  on  the  other  hand,  we  use  energy  as  the  criterion  for  keeping  zeros  (that 
is,  if  we  use  the  algorithm  as  presented  in  Table  3.1)  we  obtain  the  responses  in 
Figure  4.3  where  we  have  the  same  frequency  data  as  before,  but  now  we  use  a 
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Figure  4.1:  Concorde,  18-26  GHz,  10th  order  model,  5  and  3  sv’s  kept,  2  zeros  kept, 
0  degree  aspect  angle. 
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tenth  order  model  and  respectively  3,  5  and  7  singular  values  are  kept.  For  all  three 
cases  the  threshold  mentioned  in  Step  lb  has  been  taken  as  the  first  large  drop  in 
the  magnitude  of  the  singular  values,  which  occurs  after  the  second  singular  value 
in  all  three  cases.  It  can  be  clearly  seen  that  the  estimated  zeros  axe  very  similar 
indeed,  and  the  associated  energies  are  also  quite  similar.  If  we  compare  this  figure 
with  Figure  4.4,  in  which  we  have  dropped  the  model  order  M  from  ten  to  five,  we 
see  that  the  two  zeros  are  still  in  the  same  place,  although  the  estimated  energy  of 
the  weaker  scattering  center  differs  from  those  in  Figure  4.3.  Thus  we  see  that  the 
algorithm  presented  in  Table  3.1  is  indeed  less  sensitive  then  the  previous  algorithm 
to  the  choice  of  model  order  and  number  of  singular  values  kept. 

If  we  use  a  properly  scaled  picture  of  Concorde,  such  as  shown  in  Figure  4.5, 
to  check  for  the  geometric  relevance  of  the  two  detected  scatterers,  we  can  see  that 
the  first  peak  would  seem  to  correspond  to  the  leading  edge  of  the  wings,  and  the 
second  peak  appears  to  correspond  to  the  trailing  edge  of  the  wing  and/or  the  tail 
assembly.  Again,  the  FFT  is  provided  for  comparison  in  Figure  4.2. 

Figures  4.6  and  4.7  show  the  response  estimates  for  the  Boeing  707.  We  have 
used  N  =  30  frequencies  from  18-26  GHz  (this  corresponds  to  67-173  MHz  mea¬ 
surements  of  the  unsealed  target).  In  Figure  4.6,  the  model  order  M  =  10,  and 
the  number  of  singular  values  kept  n  is  7  and  5,  respectively.  Note  that  these  are 
the  same  parameter  choices  as  made  for  the  Concorde  in  Figure  4.3.  In  Figure  4.7, 
M  =  7  and  n  =  5  sv’s  are  kept.  In  all  these  cases  the  zero  selection  was  based  on 
the  energy  criterion,  and  5  zeros  were  kept.  Again,  we  see  that  ilic  algorithm  is 
very  insensitive  to  changes  in  the  parameters.  The  FFT  data  is  shown  in  Figure  4.8 

Comparing  the  results  with  Figure  4.9,  a  scaled  picture  of  the  Boeing  707.  we 
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Figure  4.3:  Concorde,  18-26  GHz,  10th  order  model,  3,  5  and  7  sv’s  kept,  2  zeros 
kept,  0  degree  aspect  angle. 
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Figure  4.4:  Concorde,  18-26  GHz,  5th  order  model,  3  sv’s  kept,  2  zeros  kept,  0 
degree  aspect  angle. 
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Figure  4.5:  Scaled  Picture  of  the  Concorde. 


see  that  the  middle  three  peaks  appear  to  correspond  to  the  nose,  the  antenna  and 
the  wing  complex  respectively.  The  two  outer  peaks  would  then  be  post  resonances. 
A  post  resonance  can  appear  at  a  smaller  range  than  the  main  response  because  of 
wraparound  due  to  range  aliasing. 

In  summary,  we  have  shown  that  the  algorithm  is  relatively  insensitive  to  choices 
of  parameters.  This  means  that  a  classifier  based  on  thi -  iilgorithm  could  operate 
in  a  robust  manner.  We  have  also  shown  that  the  estimates  produced  by  the  al¬ 
gorithm  seem  to  have  geometrical  relevance.  Finally,  by  comparing  the  results  for 
the  Concorde  and  the  Boeing  707  we  can  conclude  that  different  aircraft  results  in 
different  estimates. 

4.2  Tests  of  Performance  in  Noise 

The  above  comparisons  were  carried  out  using  noiseless  data.  This  section  shows 
that  the  estimates  obtained  are  relatively  insensitive  to  noise.  We  present  results 
for  the  Boeing  707  and  the  Concorde,  both  at  lOdB  signal  to  noise  ratio  (SNR).  In 
both  cases  we  have  used  a  model  order  M  =  7  and  kept  n  =  5  singular  values.  For 
the  Concorde,  we  keep  m  =  2  zeros  and  for  the  Boeing  707  we  keep  m  =  5. 

Figure  4.10  shows  the  results  of  ten  Monte  Carlo  simulations  for  the  Concorde, 
and  Figure  4.11  shows  the  results  of  the  same  experiment  for  the  Boeing  707.  The 
Figures  show  that  the  different  Monte  Carlo  simulations  lead  to  results  that  are 
similar  both  to  one  other  and  to  the  noiseless  results  (compare  with  Figures  4.3  and 
4.6  respectively).  Again,  these  are  typical  results. 

Inspection  of  the  algorithm  shows  that  this  experiment  implies  that  the  values  of 
M  and  n  can  indeed  be  chosen  independently  of  the  target,  since  the  only  difference 
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Figure  4.6:  Boeing  707,  18-26  GHz,  10th  order  model,  7  and  5  sv’s  kept, 
kept,  0  degree  aspect  angle. 
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Figure  4.7:  Boeing  707,  18-26  GHz,  7th  order  model,  5  sv’s  kept,  5  zeros  kept,  0 
degree  aspect  angle. 
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Figure  4.9:  Scaled  Picture  of  the  Boeing  707. 


between  the  parameter  values  for  the  Concorde  and  the  Boeing  707  is  the  value  of 
m,  which  does  not  influence  the  choice  of  M  or  n.  In  fact,  m  can  be  chosen  by  the 
classification  stage  after  the  signal  processing  stage  has  run  to  completion.  This 
means  that  the  proper  value  of  m  can  be  made  part  of  the  data  contained  in  the 
classifier’s  catalog. 

4.3  Other  Comparison  Tests 

This  section  describes  some  experiments  showing  the  results  of  applying  the 
algorithm  to  different  targets,  aspect  angles,  frequency  ranges,  data  lengths  and 
model  orders.  These  results  illustrate  the  range  of  applicability  of  the  algorithm. 

Figure  4.12  shows  the  Boeing  707  with  IV  =  20  data  points  taken  from  a  6.25- 
12  GHz  data  base  (this  corresponds  to  42-80  MHz  for  an  unsealed  target),  model 
order  M  —  10  and  number  of  singular  values  kept  n  —  7.  Despite  the  difference 
in  frequency  range,  the  correspondence  with  Figure  4.6  is  quite  good  (there  is  a 
difference  in  the  position  of  the  reference  point).  We  see  that  the  leftmost  peak 
in  Figure  4.6  is  missing  in  Figure  4.12;  otherwise  we  have  similar  shaped  peaks 
occurring  at  the  same  relative  spacing  in  both  Figures. 

Figure  4.13  shows  the  Boeing  707  with  N  =  40  data  points  taken  from  the 
6.25-12  GHz  range,  model  order  M  =  10,  and  n  =  7  singular  values  kept.  Except 
for  the  data  length,  these  are  the  same  choices  as  in  Figure  4.12.  Since  the  data 
length  has  been  increased  by  a  factor  of  two  while  keeping  the  bandwidth  constant, 
the  frequency  step  size  Sj  has  been  halved,  and  hence  the  maximum  unambiguous 
range  R  has  been  doubled.  As  a  result,  we  expect  to  see  the  image  of  Figure  4.12 
scaled  down  by  half  in  the  middle  of  Figure  4.13,  surrounded  by  a  flat  response 
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Figure  4.10:  Concorde,  18-26  GHz,  7th  order  model,  5  sv’s  kept,  2  zeros  kept,  lOdB 
SNR,  0  degree  aspect  angle. 
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Figure  4.11:  Boeing  707,  18-26  GHz,  7th  order  model,  5  sv’s  kept,  5  zeros  kept, 
lOdB  SNR,  0  degree  aspect  angle. 
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Figure  4.12:  Boeing  707,  6.25-12  GHz,  20  data  points,  10th  order  model,  7  sv’s 
kept,  0  degree  aspect  angle. 
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indicating  the  absence  of  target.  Comparing  Figures  4.12  and  4.13,  this  is  exactly 
what  we  see. 

In  Figure  4.14  we  see  the  Concorde,  using  N  =  20  data  points  from  the  6.25-12 
GHz  range  (48-92  MHz  for  a  full-scale  target),  model  order  m  =  10  and  n  —  7 
singular  values,  at  0  and  20  degrees  aspect  angle.  As  can  be  seen  from  the  figure, 
main  scattering  behaviour  is  similar.  The  second  large  peak  at  0  degrees  has  split 
into  two  narrower  peaks  at  20  degrees,  and  all  of  the  peaks  are  are  a  bit  closer  in 
the  20  degree  plot  ,due  to  the  apparent  shortening  of  the  target  in  the  downrange 
direction. 

Figure  4.15  shows  the  DC  10  using  N  =  20  data  points  from  6.25-12  GHz  (31- 
60  MHz  for  a  full  scale  target)  and  model  order  m  =  10  with  n  =  10  singular  values 
in  the  one  curve,  and  m  =  6  and  n  =  4  in  the  other.  Clearly,  all  major  features 
are  preserved  despite  the  changes  in  model  order  and  number  of  singular  values 
kept.  Figure  4.16  shows  a  scale  model  of  the  DC  10  which  compares  to  the  scale  in 
Figure  4.15. 

In  summary,  we  see  that  each  plane  has  its  own  characteristic  response,  which 
is  preserved  well  through  changes  in  frequency  band,  model  parameters,  and  aspect 
angle.  Moreover,  small  (20  degrees)  changes  in  aspect  angle  lead  to  predictable, 
smooth  changes  in  the  response.  As  a  result,  we  expect  to  be  able  to  identify  both 
target  type  and  target  aspect  angle  on  basis  of  the  parameters  computed  by  this 
algorithm. 
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Figure  4.13:  Boeing  707,  6.25-12  GHz,  40  data  points,  10th  order  model,  7  sv’ 


kept,  0  degree  aspect  angle. 
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Figure  4.14:  Concorde,  6.25-12  GHz,  20  data  points,  10th  order  model,  7  sv’s  kept, 
0  and  20  degrees  aspect  angle. 
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Figure  4.15:  DC  10,  6.25-12  GHz,  20  data  points,  10th  and  7th  order  model,  7  and 
4  sv’s  kept,  0  degree  aspect  angle. 
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Figure  4.16:  Scaled  picture  of  the  DC  10. 


5.  Conclusions 


We  have  developed  an  ARMA  algorithm  suitable  for  the  estimation  of  radar 
target  responses  from  frequency  domain  data.  We  have  shown  that  the  algorithm 
is  robust  with  respect  to  changes  in  the  model  parameters  and  to  changes  in  the 
frequency  band  employed  We  have  also  shown  that  the  algorithm  performs  reliably 
in  the  presence  of  noise.  Using  data  from  the  ElectroScience  Laboratory  Compact 
Range,  we  have  shown  that  different  targets  do  indeed  produce  different  response 
estimates,  and  that  varying  the  aspect  angle  of  the  target  by  a  small  amount  (20 
degrees)  leads  to  smooth  and  predictable  changes  in  the  estimated  response. 

Since  the  remaining  manual  steps  in  the  algorithm  (the  choices  of  the  number 
of  singular  values  and  zeros  to  keep)  are  easily  automated,  we  expect  that  the 
algorithm  in  its  present  form  is  suitable  to  serve  as  a  signal  processing  front  end  for 
a  radar  target  identifier.  With  such  a  configuration,  we  expect  to  be  able  to  identify 
both  the  nature  of  the  target  and  its  aspect  angle.  This  method  has  the  advantage 
that  the  signal  processing  stage  can  serve  to  both  transform  the  frequency  data 
into  geometrically  relevant  information  and  perform  data  -  ■  •mpression  (the  number 
parameters  output  by  the  algorithm  is  less  than  the  number  ol  data  points).  This 
is  important  for  the  computational  burden  imposed  by  the  classifier  stage.  The 
performance  of  such  a  two  stage  classification  procedure  needs  to  be  investigated. 

The  algorithm  itself  could  be  improved  by  studying  its  statistical  properties 
and  the  effects  of  unmodeled  dynamics  on  performance.  The  latter  is  especially 
important  in  radar  target  applications,  because  the  target  typically  contains  far 
more  scatterers  than  we  model.  In  addition  to  this,  our  model  is  an  approximation 
to  the  actual  electromagnetic  phenomena.  It  is  known  from  empirical  studies  that 
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the  singular  value  decomposition  helps  to  combat  the  effects  of  noise  and  unmodeled 
dynamics,  but  the  mechanism  behind  this  is  very  poorly  understood,  especially  for 
the  short  data  length  case,  such  as  we  have  in  the  radar  target  problem.  These 
problems  are  strong  candidates  for  further  study. 
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